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Abstract. In core-collapse supernovae, the v e and v e species may experience 
collective flavor swaps to non-electron species v x , within energy intervals limited by 
relatively sharp boundaries ( "splits" ) . These phenomena appear to depend sensitively 
upon the initial energy spectra and luminosities. We investigate the effect of generic 
variations of the fractional luminosities (7 e , 1-, l x ) with respect to the usual "energy 
cquipartition" case (1/6, 1/6, 1/6), within an early-time supernova scenario with fixed 
thermal spectra and total luminosity. We represent the constraint l e + lg + H x = 1 
in a ternary diagram, which is explored via numerical experiments (in single-angle 
approximation) over an evenly-spaced grid of points. In inverted hierarchy, single 
splits arise in most cases, but an abrupt transition to double splits is observed for a 
few points surrounding the equipartition one. In normal hierarchy, collective effects 
turn out to be unobservable at all grid points but one, where single splits occur. 
Admissible deviations from equipartition may thus induce dramatic changes in the 
shape of supernova (anti)neutrino spectra. The observed patterns are interpreted in 
terms of initial flavor polarization vectors (defining boundaries for the single/double 
split transitions), lepton number conservation, and minimization of potential energy. 
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1. Introduction 

Supernova neutrinos and antineutrinos (SN v and V) have long been studied as probes 
of both stellar and particle properties [H, [2]. In particular, core-collapse supernovae 
represent a unique laboratory for studying high-density (anti)neutrino-(anti)neutrino 
interactions and their associated flavor evolution [3]. Recent results from large-scale 
calculations [H [5] have sparked a renewed interest on this topic, which is now the focus 
of a growing literature, as reviewed in (6j[7j. An interesting outcome of all these studies 
has been the finding of neutrino flavor changes with similar ( "collective" ) features over 
extended energy ranges. 

In most supernova scenarios, the flavor evolution can often be reduced to an effective 
two-family framework (z/ e , u x ) governed by mass-mixing parameters (±Am 2 , # 13 ) [6j[7]. 
In usual notation, u x denotes either or u T , or a linear combination of z/ M and v T 
(and similarly for antineutrinos). The upper and lower sign of Am 2 refer to normal 
and inverted hierarchy, respectively. In these scenarios, collective effects are triggered 
by any nonzero value of the neutrino mixing angle #13 in inverted hierarchy (in strict 
analogy with the fall of a inverted pendulum [8]) and end up with a complete exchange 
( "swap" ) of flavor u e «-> v x [91 [H]] above a "split" energy E c dictated by v e — V e number 
conservation. 

In such effective 2u scenarios, the paradigm of a single spectral split in the neutrino 
sector — supported by detailed constructions in adiabatic approximation [11] — has been 
challenged by the observation of another, low-energy split in the antineutrino sector 
[T2| [131 EH]- Indeed, the joint occurrence of single v and V single splits has been recently 
recognized as a rather general feature in inverted hierarchy, and as a novel possibility in 
normal hierarchy [15]. Moreover, cases with double v and V splits have been identified 
as well [15]. Relaxing the assumption of energy equipartition has been instrumental in 
obtaining such multiple split cases. 

In general, the total energy luminosity L tot is distributed over six (3v + 3v) species, 

L to t = L e + Le + AL X , (1) 
or, equivalently, 

l = l e + k + 4l x , (2) 
where we have introduced the fractional luminosities 

l a = L a /L tot (a = e,e,x), (3) 

with l x = l x (we shall often ignore any distinction between v x and V x ). The usual 
assumption of "energy equipartition" amounts to take l e = 1^ = l x , namely, 

equipartition (l e , 4Z X ) = (1/6, 1/6, 4/6) . (4) 
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Figure 1. Equipartition case in inverted hierarchy. Left upper panel: flux spectra 
of v e (black, solid) and of v x (red, solid) at the end of collective effects; dotted lines 
indicate initial spectra. All spectra are formally projected at a distance d = 10 kpc. 
Left lower panel: electron flavor survival probability P ee . Right upper and lower panel: 
as before, but for antincutrinos. Vertical green lines mark the crossing energies where 
the v e (or F e ) and v x fluxes are equal. 



Figured] shows typical results at the end of collective effects (for inverted hierarchy), 
in a SN scenario endowed with the equipartition hypothesis, as also adopted in our 
previous papers (T2J [131 E] • The v e and u x differential fluxes are swapped above a split 
energy E c ~ 6 MeV. A similar swap occurs between T7 e and T7 X , but at a lower (and not 
as sharply defined) split energy E c ~ 1 MeV. No swap occurs in normal hierarchy for 
either v or V (not shown). Details of the SN scenario will be given in Sec. [3J here we just 
emphasize that significant (qualitative and quantitative) departures from such results 
can be induced by admissible deviations from equipartition, such as those considered in 
[T5] . Therefore, it is worthwhile to investigate the spectral dependence upon the initial 
SN luminosities. 

In this work, we perform an extensive phenomenological investigation in a reference 
SN scenario where all parameters are fixed, except for the fractional luminosities l a . As 
described in Sec. EJ the constraint l e + le + 4l x = 1 [Eq. (j2|)] is embedded in a "ternary 
luminosity diagram," sampled through an evenly-spaced grid of points. Details of the 
SN model used are given in Sec. [31 Numerical results for inverted hierarchy are described 
at length in Sec. HI In particular, abrupt transitions are observed from single to double 
split cases, across some regions of the ternary diagram surrounding the equipartition 
case. An interpretation of the results in terms of initial polarization vectors is offered 
in Sec. |5j Results for normal hierarchy are discussed in Sec. |6j Comments on "more 
adiabatic" scenarios are provided in Sec. [71 All our findings are briefly summarized in 
Sec. [U together with prospects for further work. 
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Figure 2. Ternary luminosity diagram. Each point corresponds to fractional 
luminosities (l e , H x ) subject to the constraint l e + I- + Al x = 1. The equipartition 
point is marked by a red triangle. The diagram is charted by an evenly-spaced grid of 
points. The small (large) shaded diamond correspond to ratios l e /l x and le/lx differing 
from 1 by less than a factor of two (four). 

2. Ternary Luminosity Diagram 

Unitarity constraints of the form hi + fi2 + h 3 = 1 can be conveniently represented 
in ternary diagrams by means of Viviani's theorem [16J, where the h^s represent the 
heights projected by any point inside an equilateral triangle of total unit height. In our 
case [Eq. (T5])], the heights are identified with the fractional neutrino energy luminosities: 
(h u h 2 , h 3 ) = (l e , 1-eAQ- 

Figure [2] shows the resulting ternary luminosity diagram. An evenly-spaced grid 
of points (at intervals 5l e = 5le = S4l x = 0.1) is superposed to guide the eye. A 
red triangle marks the equipartition point (Z e , Zg, Al x ) = (1/6, 1/6, 4/6). Departures 
from equipartition have often been considered [17]. For instance, Ref. [18] suggests to 
adopt factor-of- two- uncertainties of the kind 1/2 < l e /l x < 2 and 1/2 < lejl x < 2; 
the corresponding "allowed region" is shown as a shaded inner diamond in Fig. [2j In 
the same figure, the outer (light-shaded) diamond corresponds to a more conservative 
allowed region, embracing factor-of-four uncertainties in both l e /l x and lefl x . 
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Figure 3. Ternary luminosity diagram. Distribution of luminosities for a set of 
supernova models, as taken from the compilation in Table 7.3 of Ref. [19). 

Figure [3] shows the distribution of fractional luminosities for a set of various 
supernova simulation outputs, as taken from the compilation reported in Table 7.3 of 
Ref. [TS]. Most of the models cluster around the vertical line l-g ~ l e , except for one model 
with very high ratio h/h [20] . The overall spread in l x is rather large. It should be noted 
that the above models refer to different post-bounce times and different assumptions 
about SN (astro)physics; see [19] for a critical discussion and comparison. In any case, 
all such models are basically included in the shaded diamond which corresponds to 
a factor-of-four uncertainty in l e /l x and l^/l x , and which will be referred to as the 
"phenomeno logically interesting region." 

In this work we shall mainly focus on such restricted region for the detailed 
discussion of collective effects; however, for the sake of completeness, we have actually 
explored the full ternary diagram over all the inner points of the triangular grid, for a 
total of 36 cases (plus the equipartition one) in either hierarchies. The (computationally 
demanding) numerical investigation of all such different cases entails some assumptions 
and approximations, as described in the next Section. 
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3. Supernova Neutrino Framework 

Observations and models of SN explosions suggest a diversity of both initial and time- 
dependent features, and thus a relatively broad range of possible inputs for SN neutrino 
studies. We focus on specific SN inputs in terms of luminosities, spectra, geometry, 
frequencies, evolution equations, and selected approximations, as described below. 

3.1. Neutrino luminosities and spectra 

In previous works [321 031 EH] , inspired by the SN model in [21] , we have assumed a total 
SN binding energy Eb — 3 x 10 53 erg, and a typical luminosity decay timescale r ~ 3 s. 
At early times (t < r), where large deviations from the equipartition hypothesis may 
take place [171 [22], these choices suggests a total luminosity 

L tot ~ 10 53 erg/s , (5) 

which we fix as input, while leaving free the fractional luminosities l a = L a /L tot . 
Concerning average energies, we adopt typical early-time values [T7J, [22] , 

(E e ) = 10 MeV, (Ee) = 12 MeV, (E x ) = 15 MeV. (6) 

Finally, the initial (normalized) energy spectra Q Q (E) are assumed to be thermal, 

where £3 ~ 1.202 and j3 a is an inverse temperature parameter (/3 = 1/T) [T2], equal to 

f 0.315 MeV" 1 (a = e) , 
P a = y^- = I 0.263 MeV" 1 (a = e) , (8) 
{ a) { 0.210 MeV" 1 (a = x) , 

with 

7vr 4 

c+ = ~ 3.151 . (9) 

+ I8OC3 

3.2. Geometry 

We assume a bulb model [5] with neutrinosphere radius 

Ru = 10 km . (10) 

We also adopt the so-called single-angle approximation, where neutrino-neutrino 
interactions are averaged along a single, radial trajectory [5]. In this case, for r > R u , 
the effective v a number density per unit of volume and energy is given by [5] 

where gir) is a geometrical damping factor, 
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decreasing from unity to zero as ~ 1/r 4 . (See however [23] for a somewhat different 
function g(r) proposed in the single-angle limit.) Integration over energy provides the 
effective u a density per unit volume, 



N a (r) = / dEn a (r, E) = c{r)l a (3 a , (13) 



where 

/ \ L tot g[r) 

c(r) = ~ • (14) 

For definiteness, the absolute fluxes per unit of energy, 

Fa{E) = ' (15) 

will always be projected at a typical galactic-center distance of 

d = 10 kpc , (16) 

for both initial and intermediate spectra (at the end of collective effects). 

3.3. Frequencies 

In general, the most relevant frequencies in the context of SN neutrino collective effects 
are the vacuum oscillation frequency uj, the matter potential profile A(r), and the 
neutrino potential profile fi(r). We neglect the subdominant vacuum frequency driven 
by the smallest ("solar") 5m 2 and the tiny v il -v T matter potential difference which, as 
discussed in [2], produce only very small effects even at early times in our SN model; 
see also |24j . 

The vacuum frequency reads 

Am2 5 - 07 n -li 

w = ^iT = !Afev [km (17) 

where Am 2 = 2 x 10~ 3 eV 2 has been taken. We shall consider the energy range 
E e [0.1, 80] MeV, corresponding to u G [0.06, 50] km" 1 . Note that, as nicely shown in 
|15j . the collective v and V dynamics merge in the limit uj — > 0, namely at high E (at 
least in an effective 2u description). 

The matter potential A(r) = \/2Gi?iV e (r) (where N e is the electron number density) 
can, to a large extent, be rotated away in the equations of motion [1]. However, matter 
effects play a role, e.g., in delaying the start-up of collective effects [8], modifying some 
low-energy split features [13], generating later resonance effects for relatively large #13 
(see, e.g., [6]), and possibly inducing early decoherence [25]. By neglecting A(r), one 
trades such possible effects for a significant speed-up of the numerical calculations, which 
is a relevant gain in our systematic exploration of the ternary diagram. Therefore, in 
the spirit of Ref. |15j . we choose to neglect A from the beginning, and to focus on the 
dominant effects generated by the neutrino interaction potential /i. 
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Figure 4. Neutrino interaction potential /x(r) for our supernova model, in the radial 
range r £ [10, 10 3 ] km. The vertical spread is a consequence of slight variations of the 
curves /j,(r) over different grid points of the ternary luminosity diagram. 

The /i radial profile is given by 

/x(r) = V2G F (N + N) , (18) 

where, as we shall argue below, the appropriate v and V number densities (N and N, 
respectively) refer to an effective 2v scenario, namely 

N = N e + N x , N = N E + Nz, (19) 

so that 

H(r) = c(r)(l e p e + kfe + 2 l x (3 X ) . (20) 

Figure H] shows the corresponding family of /i(r) curves (too close to be 
distinguishable) in our adopted SN model, as obtained by varying the fractional 
luminosities l a over the 36 grid points in the ternary diagram. 

Collective effects are expected to vanish when /x(r) is somewhat below the smallest 
vacuum frequency considered, co> min ~ 0.06 km _1 . We have verified that such vanishing 
is numerically realized within r < 10 3 km at any grid point; for safety, calculations are 
extended to r = 10 3 km in all cases. 
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3.4- Flavor evolution formalism 



In the paper [14] we have performed detailed three-flavor calculations for our SN model 
at the equipartition point for different times. The results of [14] show that, up to small 
effects at the percent (or lower) level, the collective dynamics basically involve only 
two neutrino species (z/ e , v x \ while the v^^ combination orthogonal to v x remains a 
"spectator" neutrino species, and similarly for antineutrinos; see also [26]. For instance, 
a typical spectral swap can be thought to occur between v e and one of the two species 
u x , the other being unaffected during collective evolution. The spectator species enters, 
however, in the later evolution at large radius via the usual # 12 mixing effects [26J. By 
focusing our attention on the main collective effects within r < 10 3 km, we can then 
reduce the full 3z/ evolution to an effective 2v one, namely: 3z/ = (u e , v x ) © u x , where 
the relevant 2u subspace [y e , u x ) is governed by (±Am 2 , #i 3 ). As in [TJ], we choose a 
default value 

sin 2 13 = HT 6 . (21) 

We have explicitly verified, for several (non-equipartition) test cases in our SN scenario, 
that the effective 2v calculations [with /i(r) defined as in Eq. (120]) ] reproduce the results 
obtained from a full 3v approach [2] within percent accuracy. 

Concerning the 2v framework, we adopt the same notation of [12] in terms of single- 
mode Bloch vectors P and P in flavor space, appropriately extended to generic fractional 
luminosities l a . In particular, the initial conditions for the global polarization vectors 
J, J and for their difference D = J — J, are given by: 

1 ^ [dE(n e -n x )z = —- a z, (22) 



N + NJ y *' l e p e + kfe + 2l x Px 

^ dE (rig — n x ) z = e ^ x ^ x z , (23) 



N + NJ l e /3 e + kfe + 2l x p x 

D = — ^= / dE (n e -ni)i = — ^ ~ kP \ - z . (24) 

N + NJ 1 C> lJ3e + kPe + 2lxPx 

Note that the integrand in J is positive (negative) for energies below (above) the crossing 
energy where the v e and v x fluxes are equal, and similarly for J; see also Fig. [Q The D z 
component of D is a constant of motion, corresponding to u e — V e number conservation 

0- 

Finally, we remind that a potential energy V can be attached to the system, 

V ~ +z (W + W) , (25) 

where the upper (lower) sign refers to normal (inverted) hierarchy. The global vectors 
W and W are defined by energy integrals analogous to J and J in the above equations, 
but with a further factor uj in the integrand [12J. Therefore, also such vectors receive 
positive (negative) contributions below (above) the crossing energies. The role of the 
crossing energies in understanding single and multiple splits has been emphasized in 

ma. 
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Figure 5. Spectral split patterns in inverted hierarchy. Blue triangles: one "high- 
energy" v split and one "low-energy" v split. Blue squares: two v splits (and no v 
split). Yellow triangles and squares: as before, but with v and V interchanged. Red 
circles: two splits for both v and 77. 

4. Results for Inverted Hierarchy: Description 

Figure [5] shows the qualitative spectral split patterns emerging from our numerical 
exploration of the ternary luminosity diagram, in the case of inverted hierarchy. In 
the lower half of the diagram, corresponding to relatively low v x luminosity (4l x < 0.5), 
we always find one v and one V split. More precisely, the blue triangles on the right 
correspond to one dominant v split at high energy (HE) plus a minor V split at low energy 
(LE), qualitatively similar to Fig. [TJ for the yellow triangles on the left, the situation 
is reversed for v and v. As the v x luminosity increases, some \v (HE) + lu (LE) split 
cases survive (blue triangles), including the equipartition point. However, these cases 
are now flanked, on the left, by a couple of points where a double v split occurs (blue 
squares) and, on the right, by two points with a double V split (yellow squares) plus 
two points with a double split for both v and V (red circles). Therefore, a complex 
phenomenology of multiple splits [TJ] emerges in the ternary diagram, even if restricted 
to the region of phenomenological interest (the largest shaded diamond). In particular, 
small departures from luminosity equipartition can abruptly change the split patterns. 
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Figure 6. Representative grid points (numbered from 1 to 15), which cover all 
the observed split patterns. The points define a trajectory which starts close to the 
equipartition case and then spirals outwards. 

In order to better appreciate these changes, quantitative results will be shown for 
representative grid points. Figure M shows fifteen selected points, forming a trajectory 
which starts close to the equipartition case and then spirals outwards. This trajectory 
allows to explore the region of phenomenological interest, and to cover all the five 
(qualitatively different) single and double split patterns described in the comment to 
Fig. [5l The corresponding spectra are discussed next. 

Figures M and [9] show the oscillated fluxes of neutrino (left) and antineutrinos 
(right) per unit energy, at the end of collective effects (r < 10 3 km). The black and 
red solid curves refer to e-flavor and x-flavor, respectively; dotted curves represent 
unoscillated fluxes. Vertical green lines mark crossing energies, where fluxes of different 
flavor intersect, and flavor changes — if any — are unobservable. It is worth noticing that, 
in all cases, the unoscillated fluxes have qualitatively similar features, despite the large 
differences in relative luminosities: the v x flux is always lower (higher) than the v e flux, 
for energies below (above) the crossing energy E c , which ranges between 10 and 30 MeV; 
and similarly for antineutrinos, with E c in the range 5-30 MeV. Nevertheless, oscillated 
fluxes show large differences from point to point. In all panels, we indicate the point 
number (from 1 to 15) on the left, and the luminosity values (l e , l-g, 4/ x ) on the right. 




Figure 7. Neutrino (left) and antineutrino (right) energy spectra at the end of 
collective effect, for the points 1-5 in Fig. [51 The corresponding fractional luminosities 
(7 e , lx) & re indicated in each right panel. Vertical green lines mark the crossing 
energies where the v e (or v e ) and v x fluxes are equal. 



In Fig. [3, both the 1st and the 2nd point show a split pattern very similar to the 
equipartition case in Fig. [T], with a sharp HE v split and a broader, LE V split. The 
situation is suddenly reversed in the 3rd and 4th point, showing a sharp HE v split and 
a broader, LE v split. The split pattern changes again in the 5th point, where a second 
split appears for i/'s above the crossing energy, while no evident flavor change emerges 
for !7's (up to small effects in the upper tail of the spectrum). 
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Figure 8. As in Fig. but for points 6-10. 



In Fig. El the 6th point shows a split pattern similar to the last one in the previous 
figure, namely, a double v split and no V split. However, in both the 7th and 8th point, 
a double split occurs also for I7's, although the closeness of the V e and V x spectra make 
it difficult to appreciate it graphically (it will be shown in a different way below). The 
9th and 10th point show a more evident V double split, but no v flavor change (up to 
small upper-tail effects in the last case). The grid points 5-10 exhaust the cases where 
we find a double split in v and/or V spectra. 




Figure 9. As in Fig. [71 but for points 11-15. 



In Fig. [9j the cases from 11 to 13 show a pattern with a sharp HE v split and a 
broader LE V split, while the situation is suddenly reversed in the last two (14th and 
15th) cases. It should be noted that the sudden exchange of roles between v and V 
across cases 13 and 14 occurs despite the fact the unoscillated spectra and the crossing 
energies are very similar. 

In general, the results of Figs. U\M suggest the existence of a sort of "phase 
transitions" in the ternary diagram, with abrupt changes in the split patterns when 
the initial luminosities cross some "phase boundaries," which we shall try to identify in 
the next Section. 
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Figure 10. Survival probability P ee for neutrinos (left) and antineutrinos (right) in 
each of six points 5-10 of Fig. [6J where a double split is observed for neutrinos and/or 
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Finally, Fig. [TO] shows the e-flavor survival probability profile for neutrinos (left) 
and antineutrinos (right), for the six grid points numbered from 5 to 10 (see Fig. [6]), 
where double splits occur in either v or V or both (see Fig. [5]). All double splits flank 
the crossing energy (vertical green line), in agreement with the arguments of [15]. As 
also noted in [15], the values of P ee for v and V coincide at high E (i.e., for ou — > 0). 

A double-double split is evident only in cases 7 and 8. In cases 5 and 10, instead, 
there seems to be a partly developed "would-be double split" for V and u, respectively. 
In the spirit of [15], it is tempting to attribute the suppression of such emerging splits 
to the relatively high (thus unfavorable [15] ) values of the associated crossing energies; 
however, this argument is insufficient to explain the total absence of V and v double 
splits in cases 6 and 9, where the associated crossing energies are favorably lower. In 
the next Sec. El we shall develop different arguments, trying to relate (at least some 
of) such spectral features to global initial conditions and potential energy minimization. 
Additional comments will be provided in Sec. 
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Figure 11. Spectral split patterns, superposed to the three lines characterized by 
J z = (green, with J z > below the line), by J z = (blue, with J z > below the 
line), and where D z = (violet, with D z > on the right of the line). The lines mark 
transitions to different split patterns. 

5. Results for Inverted Hierarchy: Interpretation 

The authors of [J5| propose a very interesting attempt to explain analytically the 
development and the stabilization of spectral splits in terms of adiabatic invariants, 
energy minimization, and localization of the crossing energies (with details to appear in 
a forthcoming paper). Here we aim at providing another viewpoint on the interpretation 
of single and double spectral splits, with emphasis on the initial orientation of the global 
polarization vectors and on some non-adiabatic features. 

Initially, the global polarization vectors J, J and D are oriented along the z axis. 
Figure fTTI displays the curves (straight lines in our SN model) where J z = (green), 
J z = (blue) and D z = (violet) at the neutrinosphere, superposed to the same marked 
grid of points of Fig. [5j We surmise that sign changes of J z , J z and D z across these 
lines mark "phase transitions" across different split patterns in the ternary diagram, as 
argued below for the various cases reproduced from Fig. [5] (blue and yellow triangles 
and squares, red circles). 
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5.1. Cases with lu(HE) + lv(LE) split: blue triangles 

These cases, which include the equipartition one, are characterized by J z > 0, J z > 0, 
and D z > 0. In other words, both J and the smaller vector J are initially aligned 
upwards, in the unstable pendulum position [8]. As already argued in [12], the 
flavor pendulum dynamics would bring both vectors downwards, but is impeded by 
conservation of D z ; the simplest way to minimize the potential energy is then to invert 
the vector J and an equivalent fraction of the vector J, leading to a sharply defined 
"high-energy" neutrino split (whose development can be described by detailed adiabatic 
constructions [TT]). 

However, the "low-energy" antineutrino split seem to have a different, nonadiabatic 
origin, as already argued in [13] and suggested by its broader features. In the present 
work (where the effects of the matter potential A are ignored), we find that the V 
modes below the split energy E s basically experience a usual "resonant transition" 
when — uj + fiD z ~ 0, and then decouple from the collective dynamics. In the energy 
range considered (E > 0.1 MeV), this resonance is highly nonadiabatic, and leads 
to a two-level hopping probability P c — 1, corresponding to P ee ~ 1. [We have 
explicitly verified the expected decrease of Pc and P ee at lower energies, as well as 
their dependence on the dlogfi/dr profile.] In other words, the low-energy neutrino 
modes resonate with the interaction potential, and do not participate in later collective 
swaps. We find that this phenomenon roughly takes place during (and a little bit after) 
the so-called period of synchronized oscillations [27] (which is then followed by "bipolar" 
oscillations [28J). Defining fi s as the value of /i at the transition from synchronized to 
bipolar oscillations [12], the antineutrino energy modes which remain pinned upwards 
are roughly determined by the condition uj > fi s D z . In our numerical experiments, 
such condition identifies the antineutrino split energies within a factor of about 2-3; 
in particular, we find that the highest values of such energies are reached for the grid 
points closest to the D z = line. 

5.2. Cases with lu(LE) + lV(HE) split: yellow triangles 

These cases are similar to the previous one, but the roles of v and V are exchanged, 
since D z < and thus J z is larger than J z . The dominant HE split, dictated by 
-D 2 -conservation, occurs for antineutrinos. The inversion of low-energy neutrino modes 
appears to be impeded by early nonadiabatic resonances for u> > —fi s D z . 

5.3. Cases with 2u splits: blue squares 

In these two cases (identified as n. 5 and 6 in Figs. 161 and fTUl) . J is oriented upwards, 
in the unstable position, while J is oriented downwards and is thus already stable. 
However, the vector J cannot be inverted, even partially (i.e., with a single split), due 
to D z conservation: the two vectors are "locked" and remain basically equal in modulus 
and orientation at the end of collective effects, as we have explicitly verified. However, 
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one can still invert a fraction of the neutrino spectra "symmetrically" with respect to 
E c [15] , so as to preserve the modulus of J while decreasing the value of W z (and thus 
the potential energy). The double neutrino split appears thus a highly constrained way 
to lower the energy of the system, when the global polarization vectors of v and V are 
basically locked. 

Finally, we note that the point identified as n. 5 in Figs. [6] and [10] represents a 
"borderline" case, being almost on top of the line at J z = 0. This specific feature 
might be at the origin of the partial "double split," barely emerging in the antineutrino 
spectrum for this case (see Fig. [TO]) . 

5.4- Cases with 2V splits: yellow squares 

These two cases (numbered as points n. 9 and 10 in Figs. 161 and [TO]) are analogous to the 
previous ones (n. 6 and 5, respectively), but the role of v and V is exchanged. So, the 
same considerations apply, including the "borderline" features of point 10 which sits on 
top of the J z = line. 

A third case (numbered as point n. 8 in Figs. E] and [TU]) should also fall in the same 
category, being characterized by the same initial signs of the global polarization vectors. 
However, for this point (the lowest red circle in Fig. [TT]) we find double splits for both 
v and V (see Fig. [TO]) . We have been unable to explain this "exception" by sharpening 
energy minimization arguments. See, however, Sec. [7] for further comments. 

5. 5. Cases with 2u + 2V splits: red circles 

These cases include the mentioned point n. 8 (which escapes a simple explanation) 
and point n. 7 (the highest red circle in Fig. [TT]) . which can be instead more easily 
interpreted. The latter point is the only one, in our grid, having both J and J initially 
oriented downward, in the stable pendulum position. Then, the only way to further 
minimize the potential energy (without altering D z ) is to generate double splits which 
preserve the moduli of J and J while decreasing the values of W z and W z , which is 
indeed what happens at the end of the collective evolution (see Fig. [TO]) . 

5. 6. Summary for inverted hierarchy 

In inverted hierarchy, both single and double splits may appear in either the v or 
the V spectra, depending on the initial distribution of the total luminosity over the 
species (u e , V e , v x ). The split patterns observed numerically admit (with one exception) 
a relatively simple interpretation in terms of initial conditions for the global polarization 
vectors J, J and D, supplemented by arguments related to the flavor pendulum dynamics 
and its potential energy minimization. Transitions to different split patterns emerge 
when crossing the lines J z = 0, J z = and D z = in the ternary luminosity diagram. 
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Figure 12. Normal hierarchy: spectra for the single grid point where splits arc 
observed. For all other grid points, collective effects do not change the original spectra. 



6. Results for Normal Hierarchy 

In normal hierarchy, the flavor pendulum dynamics is generally trivial, since the initial 
conditions are typically very close to the stable equilibrium position [Hj. However, an 
interesting counterexample has been provided in [15] for a situation with relatively large 
v x luminosity. 

In our ternary diagram, we find that collective effects produce no flavor change in 
normal hierarchy, except for a single grid point — the one with the highest v x luminosity. 
In this peculiar case, both J and J start in the unstable position (now corresponding 
to downward orientation), and then partially reverse, in order to minimize energy while 
preserving D z . As a result, single splits emerge in both v and V spectra, as shown in 
Fig- E21 in qualitative agreement with the normal-hierarchy case discussed in [15] . 



7. Comments on a "more adiabatic" scenario 



In Sec. the proposed interpretation of the double split patterns emerging in Fig. [10] 
shows that initially stable (downward-aligned) global polarization vectors may be (cases 
7, 8) or may not be (cases 5, 6, and 9, 10) further stabilized by a double split. In view 
of the results in [15], it is worthwhile to investigate if the absence of some "would-be 
double splits" can be related to incomplete adiabaticity. To this purpose, we examine 
a hypothetical scenario, where the function n(r) is radially stretched by a factor of ten 
above the neutrinosphere, 

Kr) - A* fa + "-^r) ■ (26) 
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Figure 13. As in Fig. [101 but for the "more adiabatic" /i(r) profile in Eq. ([26]) 



The main results for such "more adiabatic" scenario are summarized in Fig. [13] 
A comparison with the analogous Fig. [TD] shows that the increased adiabaticity make 
2v + 217 splits emerge in all cases 5-10. In particular, for cases 5 and 6, we find that 
the "new" double splits provide a further (although quantitatively small) decrease of 
the potential energy for X7's, which already start from stable initial conditions (J z < 0); 
and analogously for i/'s in cases 9 and 10. In other words, the v or V systems seem to 
have "enough time" to relax further, and to minimize an already low potential energy. 
Note also that the split energies are somewhat different in Fig. [13] as compared with 
Fig. [EH especially on the low-energy side, where we have already argued (Sec. 5.1) 
that nonadiabatic resonances play a prominent role. In general, more adiabatic profiles 
for /j,(r) [such as in Eq. ( 126]) ] seem to cancel qualitative differences among the different 
double split cases commented in Sees. 5.3, 5.4 and 5.5, which merge in the single 2z/ + 2X7 
case described in Sec. 5.5. Conversely, all single-split cases remain qualitatively similar, 
as we have explicitly verified by numerical exploration of the ternary diagram. 
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Figure 14. Spectral split patterns for a scenario with increased adiabaticity (see the 
text for details). The lines, as in Fig. II H mark transitions to different split patterns. 

Figure [14] shows the resulting, "adiabatic" spectral split patterns, to be compared 
with the ones in Fig. [TTJ It appears that there are at least three distinct "phases" in the 
ternary diagram, with either single splits (at low/high or high/low energy) or double 
splits in the v and V spectra. Within the 2v + 277 region of Fig. [TH possible sub-phases 
seem to emerge as a result of incomplete adiabaticity in Fig. [TTJ 

The equipartition case is close to the crossing point of the three curves in Figs. [UJ 
and [TH Therefore, in the inverted hierarchy case, relatively small variations in the 
fractional luminosities, or in the (non) adiabatic character of the neutrino evolution, 
may induce significant variations in the observable spectral split patterns. 

Concerning the normal hierarchy case, we do not find qualitatively new features by 
using the modified function /i(r) in Eq. ( 1261) . as compared with the discussion in Sec. El 
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8. Summary and Prospects 

In core-collapse supernovae, collective effects involving high-density (anti) neutrinos 
continues to reveal surprising features. In particular, relatively small variations of the 
luminosities associated to u e , V e and v x appear to trigger abrupt changes from single- 
to double-split features in the energy spectra. 

We have investigated the effect of generic variations of the fractional luminosities 
(l e , le, 4:1 x ) with respect to the usual "energy equipartition" case (1/6, 1/6, 4/6), within 
an early-time supernova scenario with fixed thermal spectra and total luminosity. The 
constraint l e + le + 4l x = 1 has been embedded in a "ternary luminosity diagram," which 
has been numerically explored over an evenly-spaced grid of points. 

In inverted hierarchy, we have found both single- and double-split cases for either 
neutrino or antineutrino spectra, and have proposed an interpretation of these patterns 
in terms of initial orientations of the global flavor polarization vectors J and J, and 
of minimization of the potential energy, constrained by lepton number conservation 
[D z = J z — J z = const). In particular, the curves defined by J z = 0, J z = and 
D z = appear to provide "phase transition boundaries" between regions with single- 
and double-split features. The regions where at least one double split occurs (for either v 
or V) can actually merge into a single one with 2z/ + IV splits, if adiabaticity is increased 
by design. It turns out that the luminosity equipartition point is relatively close to the 
crossing point of the three separation curves; therefore, relatively small variations of 
the fractional luminosities can produce qualitatively different swaps and quantitatively 
significant changes in the final spectra, at the end of the collective favor evolution. 

On one hand, these results may provide a handle to reconstruct the original 
luminosities (Z e , le, 4l x ), if the associated spectral split patterns can be observed in future 
galactic SN explosions. On the other hand, they complicate the calculation of the diffuse 
supernova neutrino background, since past SN events may well have different relative 
luminosities and thus different spectral split features. Even within a single core-collapse 
SN event, there might be transitions from single- to double-split spectra (or viceversa) 
during the first few seconds after bounce, as a results of time-dependent changes in the 
relative luminosities, average energies, (non)adiabatic features, and possible departures 
from thermal spectra. 

The analysis of this rich phenomenology is just at the beginning, and calls for a 
deeper theoretical interpretations and for more refined numerical explorations, in order 
to get the most from observable supernova neutrino events. Concerning the results 
presented in this work, we plan to refine them by focusing on the phenomenologically 
interesting region close to the equipartition point, where further insights on the spectral 
split patterns may be expected, in both normal and inverted hierarchy, by adopting a 
more general SN framework and a denser sampling grid. 
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